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ABSTRACT 

We have obtained deep near-infrared _ftT s -band William Herschel Telescope observations of a 
sample of 15 nearby spiral galaxies having a range of Hubble types and apparent bar strengths. 
The near-infrared light distributions are converted into gravitational potentials, and the maxi- 
mum relative gravitational torques due to the bars and the spirals are estimated. We find that 
spiral strength, Q Sl and bar strength, Qb, correlate well with other measures of spiral arm and 
bar amplitudes, and that spiral and bar strengths also correlate well with each other. We also 
find a correlation between the position angle of the end of the bar and the position angle of the 
inner spiral. These correlations suggest that the bars and spirals grow together with the same 
rates and pattern speeds. We also show that the strongest bars tend to have the most open spiral 
patterns. Because open spirals imply high disk-to-halo mass ratios, bars and spirals most likely 
grow together as a combined disk instability. They stop growing for different reasons, however, 
giving the observed variation in bar-spiral morphologies. Bar growth stops because of saturation 
when most of the inner disk is in the bar, and spiral growth stops because of increased stability 
as the gas leaves and the outer disk heats up. 

Subject headings: galaxies: spiral; galaxies: photometry; galaxies: kinematics and dynamics; galaxies: 
structure 



1. Introduction 

The evolution of disk galaxies is significantly 
influenced by the main features within their disks, 
notably spirals and bars. In turn, galactic evo- 
lution also significantly influences the evolution, 
morphology, and dynamics of these features. Bars 
at the present epoch can be found in about 70% 
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of massive disk galaxies using near-infrared images 
(Sellwood & Wilkinson 1993; Knapen et al. 1999; 
Eskridge et al. 2000). The fraction is somewhat 
lower on optical images. 

The gravitational torque method (GTM, Buta 
& Block 2001) was developed to quantify bar 
strength. The idea is to transform a deprojected 
near-infrared image of a spiral galaxy into a grav- 
itational potential, and then compute the ratio of 
the tangential force to the mean radial force as 
a function of position in the plane of the galaxy. 
The potential is derived from Poisson's law af- 
ter assuming an exponential vertical density law 
and a constant mass-to-light ratio (Quillcn, Fro- 
gel, and Gonzalez 1994). The mean radial force 
represents the axisymmetric background due to 
the bulge, disk, and bar. A fully quantitative 
measure of bar strength can be defined from the 
maximum of the tangential-to-radial force ratio, as 
long ago suggested by Sanders & Tubbs (1980) and 
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Combes & Sanders (1981). The maximum force 
ratio is equivalent to the maximum gravitational 
torque per unit mass per unit square of the cir- 
cular speed. The GTM has been applied to large 
samples of galaxies in several recent studies (Block 
ct al. 2001, 2002; Laurikainen, Salo, & Rautiainen 
2002; Laurikainen & Salo 2002; Buta, Laurikainen, 
& Salo 2004, hereafter BLS; Laurikianen, Salo, & 
Buta 2004a; Laurikainen et al. 2004b). In general, 
the GTM has an advantage over other methods 
previously used to define bar strength because it 
is based on the forcing due to the bar itself and 
not just on the bar's apparent shape. The relative 
bar torque is a good measure of the importance of 
non-axisymmetric forces, which play a role in gas 
accretion, spiral arm generation, and overall disk 
evolution. 

Studies at optical wavelengths suggested that 
early-type barred galaxies are associated with 
grand design spiral structure because the bars end 
near their own corotation resonances, and at this 
point they have enough torque to drive corotating 
spiral waves outward into the disk (e.g., Elmegreen 
& Elmegreen 1985, 1989). Late-type bars could 
have different bar and spiral pattern speeds, as ob- 
served for NGC 925 (Elmegreen, Wilcots & Pisano 
1998). 

The interaction between bars and spirals is 
not well understood observationally because most 
studies have been based on optical images, which 
are confused by dust and star formation. We have 
shown, for example, that X-band spirals often 
look different from optical spirals with regard to 
pitch angle, continuity, and symmetry (Block & 
Wainscoat 1991; Block et al. 1999). Because the 
light distribution at 2.2/zm emphasizes the mass 
distribution in the old disk, the near-infrared spi- 
rals are more important for bar driving than the 
optical spirals. 

Buta, Block, & Knapen (2003, hereafter BBK) 
have shown that it is possible to use straightfor- 
ward Fourier techniques to separate the luminos- 
ity distributions of bars and spirals (see also Lind- 
blad, Lindblad, & Athanassoula 1996). The GTM 
is then used to derive maximum relative torques 
for each nonaxisymmetric feature alone, thereby 
allowing us to examine in a quantitative man- 
ner the possible correlation between bar and spi- 
ral strengths. In this paper, we apply the BBK 
method to 17 representative disk galaxies covering 



a range of early to intermediate Hubble types and 
de Vaucouleurs family classifications. Our goal is 
to investigate the relation between bar torque and 
spiral arm amplitude. 

2. Observations and Sample 

Our sample consists of the 17 galaxies listed in 
Table 1. These were selected on the basis of Hub- 
ble type (range Sab-Scd) and inclination (less than 
70°). The absolute blue magnitudes of the galax- 
ies range from —18.3 to —21.4, with an average of 
— 20.5±0.9. Thus, these are massive, high lumi- 
nosity systems, typical of the bright galaxy popu- 
lation. A range of apparent bar strengths covering 
the de Vaucouleurs families SA, SAB, and SB was 
also selected. 

Fifteen of our sample galaxies were observed 
during a total of five nights (2001 September 5, 
and 2001 October 5-8) with the INGRID cam- 
era (Packham et al. 2003) attached to the 4.2-m 
William Hcrschcl Telescope (WHT). The images 
were taken in the X-short (2.2^tm, or K s ) band, 
and have a scale of 0.241 arcsec/pix and a field 
of view approximately 4 arcminutes square. Total 
on-source exposure times average about 59 min, 
but ranged from 16 to 100 min (see Table 1). The 
observing techniques are the same as those de- 
scribed in Knapen et al. (2003). Particularly im- 
portant is the background subtraction, which was 
performed by interspersing between small blocks 
of galaxy observations several exposures of the 
same length on a blank background field. We esti- 
mate the typical sky subtraction error to be on the 
order of 0.1%-0.2% of the total background. To 
eliminate bad pixels, columns, and (in the case of 
the background frames) foreground stars, dither- 
ing was used for both the galaxy and background 
frames. We required deeper than usual exposures 
at 2/im in order to effectively measure spiral arm 
torques. The limiting surface brightness and typ- 
ical signal-to-noise ratio were estimated from our 
K s image of NGC 1530 using an independently 
measured K s flux listed in the NED. In this im- 
age, a 48 min exposure, the signal-to-noise ratio 
is 2.6 at a surface brightness level hk b = 20.0 
mag arcsec" 2 . Typical la background noise is 21.0 
mag arcsec" 2 , down to around 21.3 mag arcsec" 2 
for our longest exposures. 

In addition to the 15 WHT objects, our anal- 
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ysis includes two southern strongly-barred galax- 
ies, NGC 1365 and NGC 1433, which were part 
of previously published studies. These were ob- 
served with the CTIO Infrared Imager (CIRIM) 
attached to the 1.5-m telescope of Cerro Tololo 
Inter-American Observatory (CTIO) and have a 
scale of 1"14 pix -1 . Details of the observations of 
NGC 1365, observed at K, are provided by Regan 
& Elmegreen (1997), while those for NGC 1433, 
observed at H, are provided by Buta et al. (2001). 

3. Bar/Spiral Separation 

The GTM as used in previous studies involves 
the derivation of the maximum tangcntial-to- 
radial force ratio, a single number that charac- 
terizes the strength of nonaxisymmetric pertur- 
bations in a galaxy. BB01 inspected the two- 
dimensional ratio ( "butterfly" ) maps of tangential 
force Ft to mean radial force F r to insure that 
this maximum was mostly measuring a bar for 
their 36 galaxies. Hence, they called the force ratio 
maximum Qb- A more automated approach was 
used by Laurikainen, Salo, & Rautiainen (2002) 
and Block et al. (2002) whereby the maximum 
force ratio is derived from the radial variation of 
the maximum ratio Qt — \FT/F R\ max . Depend- 
ing on the relative importance of the bar and the 
spiral, the maximum Qt could be in the bar or 
the spiral region. For this reason, BBK proposed 
calling the maximum force ratio Q g to remove 
any ambiguity about what it represents. Since 
Q g in a barred galaxy could be affected by spiral 
arm torques, we cannot use Q g alone to assess 
whether stronger bars correlate with stronger spi- 
rals. Instead, we need to use the special technique 
outlined by BBK, a Fourier-based method that 
separates the torques due to a bar from those 
due to a spiral. However, unlike estimates of Q g , 
bar/spiral separation is not automated but de- 
pends on an iterative procedure. 

In this section, we carry out the BBK separa- 
tion procedure on each sample galaxy, and then 
compare the results with estimates based on arm- 
interarm and bar-interbar contrasts. The galaxy 
images are deprojected using mainly RC3 (de Vau- 
couleurs et al. 1991) orientation parameters, with 
revisions for three galaxies based on isophotal el- 
lipse fits (NGC 6951, 7723) or kinematic param- 
eters (NGC 1530, Regan et al. 1996). Gravita- 



tional potentials are evaluated using the method 
of Quillen, Frogcl, & Gonzalez (1994), under the 
assumptions of a constant near-IR mass-to-light 
ratio and an exponential vertical scalehcight h z 
estimated from the radial scalelength Hr and the 
type-dependent ratio hn/h z from de Grijs (1998). 
A revised lookup table for the vertical dimension 
(sec BBK) was used that did not involve any grav- 
ity softening. 

3.1. Application of the BBK Technique 

In the BBK approach, relative Fourier intensity 
amplitudes are derived from the deprojected near- 
infrared images. In the inner parts of the galaxy, 
we assume that the nonaxisymmetric amplitudes 
are due entirely to the bar, while past a radius r m 
(where m is the Fourier azimuthal index) the ob- 
served amplitudes are due to a combination of the 
bar and the spiral. We assume that the even rela- 
tive Fourier amplitudes due to the bar decline past 
the maximum at r m in the same manner as they 
rose to that maximum. A Fourier image of the 
bar is constructed and removed from the original 
image, to give a spiral plus disk image. Then the 
m=0 Fourier image is added back to the bar image 
to give the bar plus disk image. We convert each 
of these images into gravitational potentials, and 
derive maps of the ratio of the tangential force to 
the mean background radial force, the latter being 
defined by the m~0 Fourier image. From these ra- 
tio maps, we derive the maximum force ratios Qb 
due to the bar and Q s due to the spiral arms. In 
practice, we base our analysis entirely on Fourier- 
smoothed images using terms up to m=20. BBK 
fully discuss the pitfalls and uncertainties in this 
method of separation, but it is probably the most 
straightforward approach one can use for this pur- 
pose. 

Although BBK illustrated the technique using 
a fairly symmetric barred spiral (NGC 6951), our 
present sample includes one asymmetric barred 
spiral, NGC 7741, that necessitated special at- 
tention to odd Fourier amplitudes in the bar. 
NGC 7741 is a late-type (SB(s)cd) barred spiral 
with strong asymmetry in the bar, a characteristic 
which is actually fairly typical of Sd-Sm spirals (de 
Vaucouleurs & Freeman 1972). To deal with it, we 
took advantage of the fact that NGC 7741 has no 
bulge, and chose a center for the Fourier analysis 
based on the faintest discernible isophotes of the 
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bar. With such a center, the odd Fourier terms in 
the bar go to zero near the bar ends. In this cir- 
cumstance, we can define the bar by even and odd 
Fourier terms and allow for its asymmetry. We 
have not accounted for a slight asymmetry in the 
bars of two other cases, NGC 972 and NGC 7479, 
because these have bright central regions and the 
effect is less important. 

The success of bar/spiral separation depends 
on the quality of the near- infrared images. In 
our case, the WHT images arc sufficiently well- 
exposed that noise levels are minimal in the outer 
parts of the images. Very bright star-forming re- 
gions were removed from the images to minimize 
their impact. These could have a different mass- 
to-light ratio than the surrounding old disk stars, 
and could give false torque amplitudes. 

Application of the BBK method is not auto- 
matic and requires some iteration for the best 
choice of the radii r m . Plots of 7 m /7o, the am- 
plitude of the mth Fourier component relative to 
the m=0 component, were initially made for each 
galaxy and were evaluated by eye. The most reli- 
able extrapolations were chosen by examining the 
resulting bar+disk and spiral+disk images. Poor 
choices of extrapolation would sometimes leave ar- 
tifical depressions in the spiral+disk image, espe- 
cially near the ends of the apparent bar. The final 
extrapolations used are those which provided the 
cleanest-looking separations. BBK estimated the 
errors on Qb and Q s due to extrapolation uncer- 
tainties, and showed that a ±10% uncertainty in 
the choice of r2, the radius of the maximum of the 
dominant m~2 term, can lead to a ±4% uncer- 
tainty in Qb and ±10% uncertainty in Q s , at least 
for NGC 6951. We expect that this is a reasonable 
estimate of the uncertainty for the other galaxies 
in our sample as well (see section 3.2). 

We show Fourier extrapolations for the 12 sam- 
ple galaxies that required bar /spiral separation in 
Figure 8. (The ones that did not require such 
separation, NGC 908, 1058, 6643, 7217, and 7606, 
did not show a measurable bar signal.) In these 
plots, the curves are the relative intensity ampli- 
tudes, and only the lower order terms are shown 
(the actual analysis used terms up to m=20 as in 
BBK). The symbols show the mappings (interpo- 
lations and extrapolations) used for the bar. Ex- 
cept for NGC 1255, NGC 7723, and NGC 7741, 
the mappings are symmetric around a radius r m . 



For r < r m , the rising relative amplitudes are used 
as observed, while for r > r m , the amplitudes are 
extrapolated as the exact reverse of the rise. For 
NGC 1255, NGC 7723, and NGC 7741, the ob- 
served amplitudes already decline past r m , and in 
these cases a smaller amount of extrapolation was 
needed. 

These curves show some of the different char- 
acteristics of the bars in the sample galaxies. For 
example, both NGC 1365 and NGC 1433 have 
bars that can be extrapolated with relatively flat- 
topped m=2 profiles. Although the bar is nothing 
more than a weak oval in NGC 1255, it still shows 
a distinct signature in the m=2 profile that re- 
quired little extrapolation. The bar is stronger in 
NGC 7723, but it dominates the amplitudes inside 
r=25" and required little extrapolation past r m . 
In NGC 7479, the bar and strong spiral produce a 
double-humped m=2 profile, with one hump cor- 
responding to the bar and the other hump corre- 
sponding to the spiral. NGC 7741 shows an asym- 
metric m=2 profile. 

Bar/spiral separation is most reliable when the 
bar has a relatively constant position angle in the 
disk. Figure 8 shows plots of the m=2 Fourier 
phase 4>2 as a function of radius for the same 12 
galaxies as in Figure 8. The solid vertical lines 
show the radius, r(Q\,), of the bar maximum from 
the force ratio maps (see Table 2). These show 
that the assumption of a constant bar phase is 
fairly good for the more strongly barred galaxies in 
the sample, but is less good for the weakly-barred 
ones. Nevertheless, the extrapolations still pro- 
duce reasonable separations even for the weaker 
cases. 

Figure 8 shows plots of the rotation curves pre- 
dicted from the light distribution for each sample 
galaxy except NGC 7217 (see Buta et al. 1995 for 
a detailed study of this galaxy). These curves are 
derived from the m=0 component of the gravita- 
tional potential as V = \J rd<f>o jdr and are nor- 
malized to the maximum values. The curves look 
relatively typical of high luminosity spirals. 

Figures 4 and 5 show the analysis images and 
ratio maps, respectively, for the 12 sample galaxies 
requiring bar/spiral separation. Three images are 
shown for each galaxy (following BBK): the m=0- 
20 Fourier sum image, the separated bar+disk im- 
age, and the separated spiral+disk image. In all 
cases except NGC 1808, the bar or oval has been 
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rotated to the approximate horizontal position. 
For a classical barred spiral like NGC 1365, Fig- 
ure 4 shows that the extrapolations in Figure 8 
produce a fairly clean separation of the bar and 
the strong spiral. The ratio maps for NGC 1365 
in Figure 5 show a very regular pattern of alter- 
nating force maxima/minima for this galaxy. In 
the case of NGC 1433, the extrapolations in Fig- 
ure 8 also produced a clean separation, but the 
separated spiral is an oval ring-shaped pattern. 
Hence, the spiral shows a butterfly pattern sim- 
ilar to the bar. The remaining galaxies each show 
distinctive characteristics in both the images and 
ratio maps, but the main point is that the sepa- 
rations appear to be reasonable representations of 
the galaxies. In the cases of NGC 1255 and 6814, 
where the bar is nothing more than an m=2 oval, 
the four maximum points lie at about ±40°-45° 
to the bar axis, while for the more strongly-barred 
cases, these points lie much closer to the bar axis. 
As noted by BB01, higher order terms in the bar 
potential cause the maximum points to lie closer 
to the bar axis. 

Several of the more highly inclined galaxies 
in our sample, NGC 908, 972, 5033, 6643, and 
7606, required bulge/disk decompositions in or- 
der to minimize the effects of bulge "deprojection 
stretch." The idea is to model the bulge profile, 
subtract off the bulge from the 2D images, de- 
project the residual disk light, and then add the 
bulge back to the deprojected disk light. Since 
none of these five galaxies is strongly barred, one- 
dimensional decomposition techniques were ade- 
quate for our purposes. We used either the Kent 
(1986) iterative method or r< bulge and exponen- 
tial disk decompositions (e.g., Kormendy 1977). 

The bar separation is perhaps less certain for 
NGC 5033 than for the other galaxies. The bulge 
decomposition left some residual zones of lower in- 
tensity around the galaxy minor axis. A small bar 
is found in the inner regions, but its ratio map is 
affected by the decomposition uncertainties. The 
ratio map also highlights the bar- like nature of the 
spiral in this galaxy. 

Of the strongly barred spirals with bright cen- 
ters, the asymmetry in NGC 7479's bar left a resid- 
ual pattern in the spiral plus disk image. One sees 
a low intensity region on one side of the center, and 
a higher intensity region on the other side. Ignor- 
ing this asymmetry will not affect Qb too much 



in this case, since Qh is based on an average of 
the maximum points in the four quadrants, but 
will cause us to underestimate the scatter in these 
maximum points. 

Our refined procedure for NGC 7741 accounts 
for the asymmetry in the bar very well, but leaves 
a sharp edge in the area around the bar ends. This 
has very little impact on Qb, but Q s could have a 
larger uncertainty because of it. 

From the separated ratio maps, we derived 
plots of the maximum force ratio, Qt{t) = 
\F-r / Fon\ max as a function of radius r in the 
galaxy plane. The way these curves are derived 
is fully described by BBK. At each radius, we lo- 
cate the force ratio maximum in four quadrants 
and then average the results. Figure 8 shows 
the curves for 16 of our sample galaxies, includ- 
ing the ones for which no clear bar signal was 
detected. In this figure, the radii are normal- 
ized to ro (25) = D Q /2, where D Q is the Galactic 
extinction-corrected face-on blue light isophotal 
diameter at a surface brightness level of 25.0 mag 
arcsec" 2 (RC3). The results for NGC 7217 are 
not shown since little forcing due cither to a bar 
or a spiral was detected in this object. Most of 
the signal seen in the force-ratio map is likely to 
be due to bulge deprojection stretch, and has Qt 
< 0.04 everywhere. Also, NGC 7217 is the most 
bulge-dominated system in our sample (Buta et 
al. 1995), and our procedure will not account for 
this reliably since, like BB01, we have transformed 
the light distributions into potentials assuming all 
components have the same vertical scaleheight. 
For our main analysis we have set Qb ~ Qs = 
for NGC 7217. 

Table 2 summarizes the results of the separa- 
tions. Four other sample galaxies, NGC 908, 1058, 
6643, and 7606, also had little or no detectable 
bar, and we have simply set Qb=0 for these. The 
weak ovals in NGC 1255 and 6814 were easily sep- 
arated from their spirals. The remaining galaxies 
have a range of bar strengths up to <5b=0.61 in 
NGC 1530. Owing to a limited field of view, the 
spiral parameters for NGC 1808 are more uncer- 
tain than the listed errors imply. 

The radii of the maximum relative torques tend 
to be well inside the standard isophotal radius, 
for both bars and spirals. For the seven objects 
in Table 2 having Q h > 0.25, < r(Q h )/r a > = 
0.24±0.09 (s.d.) and < r(Q s )/r > = 0.42±0.09 
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(s.d.). For these same galaxies, < r(Qb)/r(Q s ) > 
= 0.57±0.13 (s.d.), so that the maxima tend to be 
well-separated on average. 

3.2. Uncertainties 

The main uncertainties in the GTM are due 
to the uncertainties in the assumed vertical scale- 
height h z , variations in the mass-to-light ratio 
(both due to dark matter or stellar population dif- 
ferences) , the adopted orientation parameters (in- 
clination, line of nodes), the bulge deprojection 
stretch, the bar thickness, the sky subtraction, 
and the galaxy asymmetry. In addition to these, 
bar/spiral separation involves uncertainties due to 
the method of extrapolation. For the more highly 
inclined galaxies in the sample, we have minimized 
bulge deprojection stretch using bulge/disk de- 
composition. 

BLS show that the uncertainties in maximum 
relative torques average about 12% due to uncer- 
tainties in the vertical scalehcight. A nonconstant 
vertical scaleheight in the bar could lead to a fur- 
ther 5% uncertainty in Qb (Laurikaincn & Salo 
2002). Uncertainties in galaxy inclinations lead 
to an inclination-dependence in the uncertainty in 
Qb or Q s . For an error of ±5° in inclination i, the 
error in maximum relative torques ranges from 4% 
for i <35° to 20% for i >60°. An uncertainty of 
±4° in line of nodes position angle <f> has less of 
an impact, ranging from 4% for i <35° to 11% 
for i >60°. Typical sky subtraction uncertainties 
on the WHT images could lead to an additional 
3.5% uncertainty in Qb and 5% uncertainty in Q s 
(BBK). Finally, uncertainties in the bar extrap- 
olations and in the symmetry assumption of the 
relative Fourier intensity amplitudes can lead to 
an additional 4% uncertainty in Qb and an 11% 
uncertainty in Q s (BBK). Table 2 lists the total 
uncertainty on each parameter from this analysis. 

For NGC 1365, the uncertainty in the orienta- 
tion parameters may be larger than we have as- 
sumed. Jorsater & van Moorsel (1995) analyzed 
an HI velocity field of NGC 1365 and concluded 
that circular motions may be reliable only be- 
tween radii of 120" and 240", where the kinematic 
position angle and inclination are relatively con- 
stant. The kinematic position angle, 220°, agrees 
well with the RC3 value of 212°. However, the 
kinematic inclination of 40° significantly disagrees 
with the nearly 58° inclination implied by the RC3 



logarithmic isophotal axis ratio of logR25—0.26. 
Since the bar of NGC 1365 is nearly along the 
galaxy's minor axis, the derived torque parame- 
ters will be sensitive to this disagreement. 

The impact of dark matter on maximum gravi- 
tational torques for the OSU bright galaxy sample 
was shown by BLS to be fairly small, especially for 
those galaxies more luminous than L»=2.5xl0 10 
Lq in the .B-band (corresponding to M^=— 20.5). 
Because the galaxies in the OSU sample are on av- 
erage luminous, massive systems, the correction to 
Q g was generally less than 10% with an average of 
about 5%. Our sample here has similar character- 
istics, with an average luminosity of Lb=L*. The 
maximum relative gravitational torques in both 
the bar and the spiral regions tend to lie in the 
bright inner regions of the galaxies, where the dark 
halo contribution is small in such luminous galax- 
ies. Thus, it is likely that dark matter has only a 
minimal impact on our results. 

The stellar mass-to- light (M/L) ratio is a sep- 
arate issue that could impact gravitational torque 
calculations. Bell & de Jong (2000) used sim- 
ple spectrophotometric evolution models of spiral 
galaxies to show that stellar M/L variations can, 
in fact, be significant, even in the if-band. They 
present simple relations between color index and 
the M/L correction for a given passband. How- 
ever, we cannot make a reliable deduction from 
their analysis of the error committed by ignoring 
these variations. Reddening can be significant in 
color index maps, and would invalidate any M/L 
corrections from the Bell & de Jong relations in 
some regions, such as bar dust lanes. Also, tests 
we have made with V — K s color index maps of 
NGC 1530 and NGC 7723 indicate that the near- 
IR spirals are decoupled from the optical spirals 
(see section 6) and arc not necessarily much bluer 
than bars, implying that the mass-to- light ratios of 
the spirals might not be very different from those 
of bars. However, interbar and interarm regions 
can be bluer than these features, implying that 
we could be underestimating Qf, and Q s slightly. 
Since our sample is defined by intermediate to late- 
type spirals, it is likely that Qf, and Q s are affected 
by stellar M/L variations in a similar manner for 
each galaxy, thus largely preserving any relation- 
ship between them. 
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4. Correlation Between Bar and Spiral 
Strength 

With our estimates of separate bar and spiral 
strengths, we can now check how well these two 
quantities correlate. Figure 7 shows Q s versus Qb 
assuming constant M/L. For Qb < 0.3, there is 
only a weak correlation with Q s , meaning that spi- 
rals as strong as Q s «0.3 are possible even in the 
absence of a bar. However, the strongest spirals 
in our sample are only associated with bars having 
Qb > 0.4. The error bars tend to be large for these 
stronger cases in part because two of the galaxies 
having Q b > 0.4 (NGC 1365 and 1530) have RC3 
inclinations close to 60° and one (NGC 7741) has 
considerable asymmetry in the bar. Figure 7 also 
suggests that among barred galaxies, stronger bars 
have stronger spirals. That is, for Qb > 0.2, Qb 
and Q s are linearly correlated. 

5. Other Indicators of Bar and Spiral 
Strength 

In addition to gravitational torques, the rel- 
ative importance of bars and spirals can be as- 
sessed in a variety of different ways. For example, 
the deprojected ellipticity of a bar is thought to 
be a good indicator of bar strength (e.g., Martin 
1995), based on the theoretical study by Athanas- 
soula (1992). Abraham & Merrifield (2000) have 
refined this idea using an automated analysis of 
two-dimensional surface brightness distributions 
and a rescaled bar ellipticity parameter called fbar 
(Whyte et al. 2002). However, bar ellipticity is 
an incomplete measure of bar strength since the 
latter also depends on the total mass of the bar 
(Laurikainen, Salo, & Rautiaincn 2002). 

Other estimates of bar importance include the 
bar-interbar contrast and the amplitude of the 
m = 2 Fourier component of the bar. Simi- 
larly, the arm strength can be assessed through 
the use of arm-interarm contrasts. We used the de- 
projected images to make these estimates for our 
sample galaxies. Arm classes (flocculent, multiple 
arm, and grand design; Elmcgreen & Elmegreen 
1982) were assigned to the galaxies according to 
their deprojected if s -band appearance, as indi- 
cated in Table 3; some assignments differ from 
their original classification in the _B-band (e.g., 
Elmegreen & Elmegreen 1987). Intensity cuts 
were made along the major axis of the discernible 



bars (using PVECTOR in IRAF) in order to de- 
termine their radial profiles. Bars with constant 
surface brightness profiles are labeled "flat" in the 
table, while bars with decreasing profiles are la- 
beled "exponential" (as discussed in Elmegreen & 
Elmegreen 1985). The lengths of the bars were 
determined using a program, SPRITE, which fits 
ellipses to isophotes via least-squares (see Buta et 
al. 1999). Bar length is taken to be the major 
axis length at maximum ellipticity, close to where 
the position angle begins to change. The max- 
imum ellipticities and bar radii r^ar (relative 
to the extinction-corrected isophotal radius from 
RC3) are given in Table 3. In addition to max- 
imum ellipticities and bar radii from the full im- 
ages, we also derived these parameters from the 
separated bar images. This is useful for those more 
strongly barred spirals where the brightest parts 
of the arms lie very close the bar ends. These are 
listed as r' bar /r and e' b in Table 3. 

Comparisons between r(Qb)/r , r(Q s )/r , and 
Tbar/fQ indicate that r(Q s ) « r^ar in several of the 
more strongly barred cases where the inner parts 
of the spirals affect the isophotal bar fit. However, 
r(Q s ) > r' bar in all cases. 

The deprojected galaxies were transformed into 
polar images, with gray scale intensities for radius 
versus azimuthal angle. The bar-interbar ampli- 
tude at Q.lrbar was estimated from intensity pro- 
files parallel and perpendicular to the bar. Also, 
azimuthal intensity cuts with a width of three pix- 
els were made from the polar images at radii cho- 
sen to be at 0.25 and 0.50 times the standard 
isophotal radius. From these profiles, the magni- 
tude differences between the bars and the interbar 
regions or between the arms and the interarm re- 
gions were determined for each galaxy, where 

Contrast(mag) = 2.5log(I pea k/ 'I disk) 

for Ip e ak — bar or arm intensity and Idisk — inter- 
bar or interarm intensity. We label the contrast 
parameters as A>, ar , A3.25, and A). 50, respectively. 
The results are compiled in Table 3. 

Fourier transforms for the m—2 components of 
the intensity cuts at 0.7rb or were made using the 
equation 

F(2) = 0£ I(e)sin(28)] 2 + [£ E( 7 W) 
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where 1(9) is the intensity at each azimuthal angle 
9 (ranging from 0° to 360°), and the sums are over 
all angles. These are also tabulated in Table 3. 

The results reveal several correlations among 
these parameters. The bar-interbar contrast A^ ar 
scales in a nearly linear fashion with Qb, but with 
a considerable scatter at large Ab ari as shown 
in Figure 8a. The bar ellipticity also correlates 
with bar torque, as shown in Figure 8b, although 
the relation is not linear. This is fully consistent 
with the results of Laurikaincn, Salo, & Rauti- 
ainen (2002), who derived this same correlation 
from images of more than 40 galaxies from the 
Two Micron All-Sky Survey (2MASS). For Q b < 
0.4, the correlation between ej, and Qb is fairly 
good, while for Qb > 0.4, the correlation is weaker 
in the sense that small changes in can corre- 
spond to large changes in Qb- The correlation 
is only slightly improved when e' h is used. The 
Fourier strength parameter F(2) correlates with 
Qb (with some scatter at large torque), as shown 
in Figure 8c. These results indicate that the bar- 
interbar contrast, ellipticity, Qb, and m=2 Fourier 
component all provide reasonable measures of bar 
strength, especially for the weaker bars. 

The bar-interbar contrast increases with in- 
creasing bar radius, as shown in Figure 8d, where 
Abar is plotted versus the fractional bar radius 
rb ar /rQ. Furthermore, the bar torque is stronger 
for fiat bars (average value of Qb is 0.43±0.14 for 6 
objects) than for exponential bars (0.28±0.16 for 
4 objects). These results are consistent with pre- 
vious K-b&nd studies, which indicated that flat 
bars are both longer and stronger than exponen- 
tial bars (Elmcgreen & Elmegreen 1985, Regan 
& Elmegreen 1997). The previous studies also 
found that earlier spiral types tend to have flat 
strong bars and later types tend to have expo- 
nential weak bars, although a correlation with 
Hubble type is not obvious in the present small 
sample. However, correlations of maximum rela- 
tive gravitational torques with Hubble type have 
been established by Laurikainen, Salo, & Rauti- 
ainen (2002), BLS, and Laurikainen, Salo, & Buta 
(2004). Early-type bars and spirals tend to be 
weaker than late-type bars and spirals owing to 
the dilution effect of bulges in early- types. In 
the present sample, the strongest bar is that in 
NGC 7741, which also lacks any significant bulge. 
Thus, the scatter in the correlation between bar 



strength and bar length could be influenced by 
the importance of the bulge. 

The arm-interarm contrast A0.25 scales with 
the arm-interarm contrast A0.50 (not shown) but 
is slightly weaker, which is consistent with pre- 
vious findings that the arm-interarm contrast in- 
creases with radius to about mid-disk (Elmcgreen 
& Elmegreen 1985). There is a weak correlation 
(not shown) between arm strength and Arm Class, 
with grand design galaxies having stronger spirals 
than flocculent galaxies. The arm-interarm con- 
trast increases linearly with the maximum relative 
spiral torque Q s , as shown in Figure 8e. The cor- 
relation is poorer for ^0.50- 

6. Dust-Penetrated Classification 

We have applied the dust-penetrated classifi- 
cation scheme of Block & Puerari (1999) to each 
of our sample galaxies. On the basis of depro- 
jected near-infrared images, evolved stellar disks 
may be grouped into three principal dust pene- 
trated archetypes: those with tightly wound stel- 
lar arms characterized by pitch angles at K s of 
~ 10° (the a class), an intermediate group with 
pitch angles of ~ 25° (the (3 class) and thirdly, 
those with open spirals demarcated by pitch an- 
gles at K s of - 40° (the 7 bin). To take full 
cognizance of the duality of spiral structure and 
decouplings between gaseous and stellar disks, it 
has been demonstrated (e.g. Block & Wainscoat 
1991) that we require two classification schemes - 
one for the Population I disk, and a separate one 
for the Population II disk. A near-infrared classi- 
fication scheme can never replace an optical one, 
and vice- versa, because the current distribution of 
old stars strongly affects the current distribution 
of gas in the Population I disk. 

For dust-penetrated classification, logarithmic 
spirals of the form r = r Q cxp(— m^/p max ) (see 
Danver 1942) arc employed in the decomposition. 
The amplitude of each Fourier component is given 
by (Schroder et al., 1994) 

_ S^S^jOn r, 9) exp [~i(m9 +p In r)) 



where r and 9 are polar coordinates, /(In r, 9) is 
the intensity at position (In r, 9) , m represents the 
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number of arms or modes, and p is the variable as- 
sociated with the pitch angle P of a given m-mode, 

7YI 

defined by tanP = . We select the pitch 

Pmax 

angle of the dominant m-mode to define the dust- 
penetrated pitch angle classes described above. Of 
course, the spiral arms in barred galaxies often de- 
part from a logarithmic shape. As noted by Block 
et al. (2001), the arms may break at a large an- 
gle to the bar and then wind back to the other 
side, as in a "pseudoring." We minimize the im- 
pact of non-constant pitch angle due to rings or 
pseudorings angles by excluding from our analysis 
the bar regions of the galaxies in question, and by 
restricting the fits to a limited range in radius. 

To illustrate what the Fourier method is ex- 
tracting, we show in Figure 10 m=2 inverse 
Fourier transform contours superposed on our de- 
projected K s image of NGC 1530. This galaxy 
has a very open two-armed spiral breaking from 
near the ends of its bar. The dominant m-mode is 
in this region, and the contours are a reasonable 
representation of the pitch angle of the arms. 

Table 4 summarizes the results of our pitch an- 
gle analysis. The radius ranges used for the fits 
are listed in column 5, except for NGC 1365 and 
1433 whose classifications are from Buta & Block 
(2001). Wc continue to find a ubiquity of low or- 
der m = 1 and m = 2 modes in this sample, con- 
sistent with earlier studies (e.g. Block & Puerari 
1999). The most uncertain classifications are for 
NGC 1808 and NGC 5033, where the field of view 
limits the part of the main spiral that we can see. 

Under very special circumstances, dominant 
modes with m greater than 2 may develop within 
the modal theory of spiral structure. Block et 
al. (1994) and Bcrtin (1996) hinted that in a 
gas-rich system, some dominant higher-m modes 
should develop, and this might also induce some 
response in the stellar disk, for the stronger cases. 
Furthermore, non-linear modes may couple and 
again give rise to higher-m structures: m = 2 
and m = 1 combine to give m — 3, m — 2 and 
m = 2 combined to give m = 4, etc. (Elmcgreen, 
Elmegreen, & Montenegro 1992; Block et al., 1994 
and G. Bertin, private communication). We use 
the terminology H3 and H4 for these third and 
fourth harmonics, to assist with easy visualization 
of their evolved disk morphologies. 

Finding examples of true three- and four-armed 



spiral galaxies in the near-infrared is a great chal- 
lenge, requiring the investigation of many galaxies 
to encounter one or two unambiguous cases. A 
good example of an m = 4 stellar disk in the 
present study is NGC 6814, while a previously 
published example is ESO 566—24 (Buta et al. 
1998; Rautiainen et al. 2004). 

We also note that early Hubble type galaxies 
(e.g. NGC 972, Sab) can present very wide open 
arms in the near-infrared; NGC 972 belongs to the 
7 class. Three other examples are NGC 1808 (of 
RC3 type Sa but of dust penetrated arm class 7), 
NGC 1530 and NGC 1365 (both Hubble type SBb; 
near-infrared arm class 7). Conversely, it is also 
possible for spirals, classified as late type in the 
optical regime, to have arms in the near-infrared 
which are not very wide open (e.g. NGC 908 and 
NGC 1058 are of Hubble type Sc, but both belong 
to the more tightly wound bin: see Table 4). 

The present study continues to illustrate de- 
couplings between gaseous and stellar disks and 
the great advantage of dust-penetration. In the 
optical, NGC 972 is classified as flocculent by 
Elmegreen & Elmegreen (1987; see also Table 3). 
The photograph reproduced by Sandage & Bedke 
(1994) in panel 148 fully supports the flocculent 
designation. Sandage & Bedke (1994) note that 
the optical image of NGC 972 "is dominated by the 
heavy dust lanes crossing the near-side of the high- 
surface-brightness bulge". In the near-infrared, 
NGC 972 presents a principally two-armed spi- 
ral; the dominant harmonic is m=2. Its grand 
design, two-armed, evolved stellar disk appears to 
be completely decoupled from its fleece-like floc- 
culent Population I gas disk. While many opti- 
cally flocculent spirals may still present a floccu- 
lent appearance in the near-infrared (Elmegreen 
et al. 1999), decouplings of the two components, 
when present, are indeed very striking. 

NGC 1808 is famous for its dust lanes which 
appear to radiate almost perpendicularly to the 
major axis. Sandage & Bedke (1994) note in their 
panel 193, that "the central region of NGC 1808 
provides what appears to be direct evidence of a 
galactic fountain composed of narrow dust lanes 
perpendicular to the plane." The near-infrared 
imaging successfully penetrates the dust in this 
SABa spiral, yielding an inner pair of grand de- 
sign, wide-open arms (the dominant harmonic in 
NGC 1808 is m=2). In this case, bar/spiral sep- 
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aration has the added uncertainty that the bar 
position angle is not constant (see Figure 8) and 
treats the wide-open spiral as part of the bar. Also 
in this case, there is a strong coupling between the 
gas and star-dominated components, although the 
inner near-infrared morphology is far more regu- 
lar. 

We use Table 1 of BB01 to define three gravita- 
tional torque classes: the total torque class based 
on <3 g , the bar torque class based on Qb, and 
the spiral torque class based on <5 S . The dust- 
penetrated type (DP-type), written in the form 
"harmonic class-pitch angle class-bar torque class" 
following BB01, is given in the last column of Ta- 
ble 4. Except for NGC 7741, the highest spiral 
torque classes all belong to the highest pitch angle 
class, 7. Six of the galaxies could not be assigned 
a DP-type. 

7. Discussion 

The theory of bar-driven spirals is some- 
what controversial. Kormendy & Norman (1979) 
showed that of 33 galaxies having differential ro- 
tation and global spiral structure, nearly 80% are 
barred and likely to have their spirals driven by the 
bars. The implication was that bars and spirals 
would have the same pattern speeds. However, 
Sellwood & Sparke (1988) showed that in n-body 
numerical simulations where an initially unstable 
disk forms both a bar and a spiral, the spiral has 
a lower pattern speed than the bar. Sellwood & 
Wilkinson (1993) suggested that the quadrupolc 
moment of a realistic bar falls off too rapidly to 
drive the spirals well beyond the ends of the bar. 
In this circumstance, the spiral could be an in- 
dependent instability, or a driven response to a 
resonance interaction between the bar and spi- 
ral. Tagger et al. (1987) and Yuan & Kuo (1997, 
1998) have shown that resonance interactions be- 
tween the bar and the spiral can be an efficient 
mechanism of wave amplification. 

The question of whether spirals and bars coro- 
tate can be answered in some cases by the ob- 
served alignment of resonance rings. Inner rings 
and pseudorings in SB galaxies tend to be aligned 
parallel to bars, while outer rings and pseudor- 
ings may be aligned parallel or perpendicular to 
bars (Kormendy 1979; Schwarz 1984; Buta 1986, 
1995). The morphologies of the various ring types 



suggest they are driven by bars: inner rings of- 
ten have pointy oval shapes, indicating a con- 
nection with the inner 4:1 ultraharmonic reso- 
nance (UHR, Buta & Combes 1996; Salo et al. 
1999), while outer pseudorings have three distinc- 
tive morphologies that suggest a connection to 
the outer Lindblad resonance (OLR). These are 
the R' 1; R2, and RiRJ, morphologies known as the 
"OLR subclasses" (Buta 1985, 1986, 1995; Buta 
& Crocker 1991). The first two types were pre- 
dicted by Schwarz (1981) from test-particle simu- 
lations as being possibly due to differences in the 
gas distribution around the OLR, while Byrd et al. 
(1994) showed that the three morphologies may be 
connected in an evolutionary manner. 

Since resonance rings probably develop from 
the secular evolution of spiral patterns (Schwarz 
1981; Byrd et al. 1994; Rautiainen & Salo 2000), 
the frequent alignments suggest that the spirals 
that form the rings are driven by the bars at the 
same pattern speeds. Buta & Combes (1996) fur- 
ther argue that the co-existence of several rings in 
the same galaxy, showing shapes and alignments 
that are compatible with periodic orbits near res- 
onances (e.g., NGC 3081, Buta & Purcell 1998), 
speaks against the existence of several patterns 
with different pattern speeds. Salo et al. (1999) 
show that a model that assumes the bar and spiral 
co-rotate can fit the morphology and the velocity 
field of the early- type ringed galaxy IC 4214. 

The situation is a little different with the cen- 
tral regions of barred galaxies. Here it is possible 
for the primary bar to induce enough mass flow 
into the central kiloparsec to cause an inner bar 
instability, leading to formation of a secondary bar 
with a different pattern speed (Pfcnnigcr & Nor- 
man 1990). Support for the idea of a different pat- 
tern speed comes from studies of relative primary 
bar, secondary bar position angles (e.g., Buta & 
Crocker 1993; Wozniak et al. 1995). 

Rautiainen & Salo (2000, hereafter RS) have re- 
considered ring formation in n-body models with 
dissipatively colliding test particles. These mod- 
els show the same features as previous test-particle 
models except that with the self-gravitating stel- 
lar disk, additional spiral modes can develop in 
the same manner as in Sellwood & Sparke (1988). 
RS consider how these modes, which often have a 
pattern speed less than that of the bar, affect the 
morhology and formation of outer rings and pseu- 
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dorings. The OLR subclasses of outer rings and 
pseudorings still form in these models, but cyclic 
changes between the morphologies can occur due 
to the influence of the modes with a lower pattern 
speed. If the influence of the lower pattern speed 
modes is significant enough, RS show that mis- 
alignments may be possible. Observed misalign- 
ments between inner and outer rings and bars are 
rare but do occur (Buta 1995; Buta, Purcell, & 
Crocker 1995). 

The only prominent ringed galaxies in our sam- 
ple are NGC 1433 and NGC 7723. Of these, NGC 
1433 shows many features in common with test- 
particle simulations, including a highly elongated 
and aligned inner ring likely connected with the in- 
ner UHR, an R^ outer pseudoring likely connected 
with the OLR, a misaligned nuclear ring, and two 
secondary spiral arcs, or plumes, in the vicinity 
of the outer UHR (Buta et al. 2001). Elmegreen 
& Elmegreen (1985) argue that a bar like that in 
NGC 1433 could still be growing through the slow 
loss of angular momentum to the stellar spiral. 
In NGC 7723, deprojection of a blue-light image 
indicates that the inner ring is a nearly circular 
feature from which a prominent multi-armed spi- 
ral pattern emerges. This galaxy does not fit into 
available test-particle simulations nearly as well as 
does NGC 1433. 

We can see in the plots in Figure 8 that some of 
the bar torques in our sample are weak in the spi- 
ral region, as predicted by Sellwood & Wilkinson 
(1993). This is the case for NGC 972, 1255, 5033, 
and 6814. This may also be true in the Milky Way. 
Our Galaxy has a bar about 1.5 kpc long which 
generates the 3 kpc arm and perhaps makes the 
bulge, which is clearly bar-like and about the right 
length (see, e.g, Cole & Weinberg 2002). Corota- 
tion is at «3 kpc. The spirals outside this radius, 
which include the local spiral, the Sgr-Carina spi- 
ral, and the Perseus spiral, may be independent 
of the bar. This means their pattern speed is un- 
related to the bar pattern or perhaps related in 
a complex way through resonant excitation. In 
either case, the bar strength and the spiral arm 
strength should not be simply related. For res- 
onance excitation, even a weak perturbation can 
cause a strong response. 

Another test of the corotation of bars and spi- 
rals concerns the place in azimuth where the inner 
part of the spiral lies relative to the bar. If the pat- 



tern speeds are different or the bar and spiral are 
unrelated, then the distribution of azimuths for 
these inner-limit spiral points should be uniform. 
If all spirals end in their inner regions at the same 
angular distance from the end of the bar, then the 
bars and spirals are probably related. To provide 
a preliminary answer to this question, we exam- 
ined the twelve barred galaxies in our sample and 
visually estimated 9s, the angle between the bar 
and the inner-limit spiral points, for each promi- 
nent spiral arm in each galaxy. Table 3 compiles 

< \8s\ >, the average over the prominent arms, 
and Figure 9 shows how this average correlates 
with Qb- The plot shows that for all but two of 
the galaxies, the spirals appear to begin within 20° 
of the bar axis. The two objects having < \8s\ > 
> 40° arc NGC 5033 and 7723, the latter having 
three spiral arms forming part of the inner ring. 
Thus, for none of the most strongly barred galax- 
ies in our sample does the near-IR spiral appear 
to begin at a large intermediate angle to the bar 
axis. However, in blue light, the situation can be 
different. For example, both NGC 1530 and 6951 
could be interpreted as having serious mismatches 
between the ends of the bar and the beginnings 
of the spiral. On the other hand, the weak inner 
pseudorings in these two objects are made from 
arms that break from the ends of the bar and wrap 
around the other ends. 

8. Conclusions 

From a sample of 17 bright spiral galaxies hav- 
ing a range of bar and spiral morphologies and 
Hubble types, we find the following from i\ s -band 
images: 

1. Bars and spirals can be effectively separated 
using a Fourier-based image analysis technique 
(BBK). From this separation, we can quantify the 
bar and spiral strengths in terms of tangential-to- 
radial force ratios. In the sample, bar strengths 
Qb range from to 0.75 while spiral strengths Q s 
range from to 0.46. 

2. The spiral strength Q s correlates with Qb m 
a nonlinear fashion: spirals form in non-barred 
galaxies so Q s is independent of Qb when Qb 

< 0.3, but spiral strengths increase linearly with 
bar strengths when bars are present, for which 
Qb >^ 0.3. The effect is not an artifact of stellar 
M/L variations or dark matter in these high lu- 
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minosity spirals. If bars and spirals grow together 
in a global disk instability, then the linear relation 
between their strengths imply they have about the 
same growth rates. 

3. The bar-interbar contrast Af, ar and the depro- 
jected bar ellipticity correlate with Qb, with 
a smaller scatter for e&. The arm-interarm con- 
trast also correlates with Q s , with a better cor- 
relation for the contrast measured at a radius of 
0.25r25 compared to 0.5r25- Strong bars tend to 
be longer relative to the galaxy isophotal radius. 
This length-strength correlation may be the result 
of bar elongation over time as more and more stel- 
lar orbits join the bar potential (see simulations in 
Combes & Elmegreen 1993). 

4. Dust-penetrated classifications indicate that 
the strongest bars with the strongest spirals tend 
to have the most open spiral arms. Their pitch an- 
gles are 40° or more and their pitch angle class is 7. 
The openness of a spiral arm depends on the rel- 
ative size of the Toomre (1964) length, 2ttGT,/k 2 
for total disk column density £ and epicyclic fre- 
quency k. This length is approximately the sep- 
aration between the arms. If this length is large 
compared to the radius, then the arms are very 
open. The ratio of the Toomre length to the 
galaxy radius is approximately the ratio of the 
disk mass to the halo mass. Thus open spirals 
have relatively massive disks (see also Bcrtin ct 
al. 1977; Elmegreen & Elmegreen 1990). Our ob- 
servations imply that the most massive disks have 
the strongest bars and spirals. This suggests that 
bars and spirals form together in a global disk in- 
stability. 

5. The correlations between the maximum rela- 
tive gravitational torques of the bars and spirals 
and between a number of other measures of bar 
and spiral amplitudes, in addition to the small an- 
gular displacements between the ends of the bars 
and the inner parts of the spirals, and the align- 
ments of most resonance rings with their bars, all 
suggest that the bars and spirals in most of our 
sample have shared the same pattern speeds for 
cosmological times. In these cases, the bar and 
spiral corotation radii are the same. Other cases 
with more irregular ring structures or no evident 
correlations between bar and spiral strengths or 
alignments could have their spirals excited either 
independently of the bar, or excited at a higher 
order resonance, giving the bars and spirals differ- 



ent pattern speeds and no correlation in structure. 
The Milky Way is apparently in this latter cate- 
gory. 

These conclusions do not necessarily apply to 
SBO galaxies, which have strong bars and weak 
(or no) spirals. This difference from our present 
sample illustrates the changing morphology of bar- 
spiral patterns over time. Our results here suggest 
that the growth phase of a strong bar-spiral pat- 
tern is the result of a combined instability having 
one pattern speed and one growth rate in a relative 
massive disk. This growth may be spontaneous or 
it may follow an encounter with another galaxy. 
After some time, the bar-spiral growth should slow 
or stop, but it does this for different reasons in the 
bar and spiral regions. Bars stop growing when 
they saturate to very large strengths, placing most 
of the disk stars within the inner 4:1 resonance in- 
side the bar. The inner stellar population is very 
hot when this happens, because the stellar orbits 
are highly non-circular. Spirals stop growing when 
the stellar disk near corotation also heats up, but 
for the spirals, the high stellar velocities are more 
random than for the bar, removing the spiral pat- 
tern. The high stellar dispersion also requires a 
nearly complete conversion of gas into stars so that 
gas dissipation is absent and young stars with low 
dispersions no longer form. The result is a rela- 
tively strong, squared-off bar with little dust and 
star formation structure in the disk and very weak 
or no spiral arms, making the SBO class. 
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Table 1 

Properties of the Sample Galaxies' 1 



Galaxy 


Type 


logD„ 


logi?25 




A 

^ivipc ) 


M° 


source 


Exp. Time 
(mill) 


1 


2 


3 


4 


5 


6 


7 


8 


9 


NGC 908 


SA(s)c 


1.78 


0.36 


10.35 


17.8 


-20.9 


WHT-INGRID 


73 


NGC 972 


Sab 


1.55 


0.29 


11.48 


21.4 


-20.2 


WHT-INGRID 


64 


NGC 1058 


SA(rs)c 


1.51 


0.03 


11.55 


9.1 


-18.3 


WHT-INGRID 


62 


NGC 1255 


SAB(rs)bc 


1.62 


0.20 


11.26 


19.9 


-20.2 


WHT-INGRID 


60 


NGC 1365 


SB(s)b 
(R')SB(r)ab 


2.05 


0.26 


9.93 


16.9 


-21.2 


CTIO-CIRIM 


40 


NGC 1433 


1.81 


0.04 


10.64 


11.6 


-19.7 


CTIO-CIRIM 


30 


NGC 1530 


SB(rs)b 


1.72 


0.28 


11.42 


36.6 


-21.4 


WHT-INGRID 


48 


NGC 1808 


(R)SAB(s)a 


1.81 


0.22 


10.43 


10.8 


-19.7 


WHT-INGRID 


56 


NGC 5033 


SA(s)c 


2.03 


0.33 


10.21 


18.7 


-21.2 


WHT-INGRID 


16 


NGC 6643 


SA(rs)c 


1.60 


0.30 


11.14 


25.5 


-20.9 


WHT-INGRID 


63 


NGC 6814 


SAB(rs)bc 


1.54 


0.03 


11.32 


22.8 


-20.5 


WHT-INGRID 


60 


NGC 6951 


SAB(rs)bc 


1.68 


0.08 


10.71 


24.1 


-21.2 


WHT-INGRID 


60 


NGC 7217 


(R)SA(r)ab 


1.63 


0.08 


10.53 


16.0 


-20.5 


WHT-INGRID 


64 


NGC 7479 


SB(s)c 


1.63 


0.12 


11.22 


32.4 


-21.3 


WHT-INGRID 


59 


NGC 7606 


SA(s)b 


1.73 


0.40 


10.88 


28.9 


-21.4 


WHT-INGRID 


58 


NGC 7723 


SB(r)b 


1.55 


0.17 


11.57 


23.7 


-20.3 


WHT-INGRID 


48 


NGC 7741 


SB(s)cd 


1.65 


0.17 


11.43 


12.3 


-19.0 


WHT-INGRID 


100 



a Col. 1: galaxy name; col. 2: dc Vaucouleurs revised Hubble type (RC3); col. 3: log of corrected isophotal 
diamatcr (units of O'l) at (j, B — 25.0 mag arcsec -2 (RC3); col. 4: log of iosphotal axis ratio at fi B — 25.0 
mag arcsec -2 (RC3); col. 5: total corrected blue light apparent magnitude (RC3); col. 6: distance in Mpc 
(Tully 1988); col. 7: absolute blue light magnitude: col. 8: source of image (WHT— William Hcrschel Telescope; 
CTIO— Ccrro Tololo Inter-American Observatory); col. 9: total on-source exposure time in minutes 



Table 2 
Derived Torque Parameters' 1 



Galaxy 


h z 


Qb 


Q a 


Qg 


r(Q b ) 


r(Q 8 ) 


r(Qg) 


r(Q b )/r 


r(Q s )/r 


r(Q g )/r„ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


NGC 908 


379 


(0.00) 


0.28±0.10 


0.28 




73.0 


73.0 




0.40 


0.40 


NGC 972 


372 


0.22±0.06 


0.21±0.06 


0.24 


9.0 


29.0 


8.0 


0.08 


0.27 


0.08 


NGC1058 


100 


(0.00) 


0.09±0.02 


0.09 




18.0 


18.0 




0.19 


0.19 


NGC1255 


482 


0.07±0.01 


0.15±0.04 


0.17 


4.5 


23.5 


76.0 


0.04 


0.19 


0.61 


NGC1365 


1587 


0.40±0.11 


0.36±0.11 


0.49 


118.5 


178.0 


164.0 


0.35 


0.53 


0.49 


NGC1433 


581 


0.37±0.06 


0.23±0.05 


0.43 


68.0 


106.5 


69.0 


0.35 


0.55 


0.36 


NGC1530 


777 


0.61±0.16 


0.42±0.16 


0.73 


45.0 


65.5 


49.0 


0.29 


0.42 


0.31 


NGC1808 


600 


0.22±0.04 


0.10±0.02 


0.24 


70.5 


118.0: 


73.5 


0.36 


0.61: 


0.38 


NGC5033 


841 


0.07±0.02 


0.12±0.03 


0.12 


8.5 


29.5 


29.5 


0.03 


0.09 


0.09 


NGC6643 


383 


(0.00) 


0.21±0.08 


0.21 




28.0 


28.0 




0.23 


0.23 


NGC6814 


375 


0.07±0.01 


0.09±0.02 


0.10 


10.0 


44.5 


15.0 


0.10 


0.43 


0.14 


NGC6951 


640 


0.28±0.04 


0.21±0.06 


0.34 


31.5 


57.5 


43.0 


0.22 


0.40 


0.30 


NGC7217 


566 


(0.00) 


(0.00) 


<0.04 














NGC7479 


494 


0.59±0.10 


0.46±0.12 


0.71 


31.0 


50.0 


45.0 


0.24 


0.39 


0.35 


NGC7606 


898 


(0.00) 


0.08±0.03 


0.08 




42.5 


42.5 




0.26 


0.26 


NGC7723 


508 


0.30±0.05 


0.12±0.02 


0.31 


16.0 


33.0 


16.0 


0.15 


0.31 


0.15 


NGC7741 


274 


0.74±0.22 


0.35±0.07 


0.77 


15.0 


46.0 


15.0 


0.11 


0.34 


0.11 



a Col. 1: galaxy name; col. 2: vertical scalchcight in pc; col. 3: relative bar torque parameter (bar strength); col. 4: relative 
spiral torque parameter (spiral strength) ; col. 5: gravitational torque parameter (total nonaxisymmctric strength) ; col. 6: radius 
(arcsec) of maximum average relative bar torque; col. 7: radius (arcsec) of maximum average relative spiral torque; col. 8: radius 
of maximum average relative total torque; cols. 9-11: same radii as in cols. 6-8, relative to r — D /2, the extinction-corrected 
de Vaucouleurs isophotal radius. 
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Table 3 

Spiral and Bar Contrasts, Dimensions, and Shapes' 1 



Galaxy 
1 


Arm Class 
2 


^0.25 

3 


^0.50 

4 


A bar 
5 


Tbar/ro 

6 


7 


r 'bar/ r 
8 


4 
9 


F(2) 
10 


< |0 S | > 
11 


NGC 908 


M 


0.59 


0.79 
















NGC 972 


F 


0.96 


1.27 


0.53 


0.11 


0.51 


0.10 


0.51 




3:±1: 


NGC 1058 


F 


0.38 


0.52 
















NGC 1255 


M 


0.56 


0.67 


0.20 


0.06 


0.18(exp) 


0.06 


0.18 




18±17 


NGC 1365 


G 




2.10 


2.05 


0.53 


0.76(flat) 


0.50 


0.71 


0.40 


9±4 


NGC 1433 


G 




1.78 


1.64 


0.54 


0.70(flat) 


0.37 


0.65 


0.33 


8±2 


NGC 1530 


G 


1.67 


3.01 


2.35 


0.43 


0.74(flat) 


0.33 


0.67 


0.38 


15±5 


NGC 1808 


G 


0.39 


0.96 


0.69 


0.37 


0.52(exp?) 


0.37 


0.52 


0.06 




NGC 5033 


M 


0.65 
















67±5 


NGC 6643 


M 


0.73 


0.67 
















NGC 6814 


M 


0.38 


0.76 


0.37 


0.12 


0.25(cxp) 


0.16 


0.27 




4±3 


NGC 6951 


G 


1.25 


1.63 


1.15 


0.33 


0.64(flat) 


0.25 


0.58 


0.22 


12±7 


NGC 7217 


F 


0.19 


0.34 
















NGC 7479 


G 


2.03 


1.46 


1.47 


0.39 


0.71(flat) 


0.23 


0.68 


0.93 


7±7 


NGC 7606 


G 


0.43 


0.33 
















NGC 7723 


M 


0.54 


0.57 


1.02 


0.18 


0.61(flat) 


0.18 


0.61 


0.11 


51±22 


NGC 7741 


F 


1.35 


0.94 


1.34 


0.14 


0.73(cxp) 


0.14 


0.73 


0.32 


7±4 



a Col. 1: galaxy name; col. 2: spiral Arm Class (Elmcgrccn & Elmcgrccn 1987, where F— flocculent, M— multi-armed, 
G— grand-design); cols. 3,4: the near-IR arm-interarm contrast at 0.25 and 0.50 times the RC3 isophotal radius; col 5: 
the near-IR. bar-intcrbar contrast at 70% of the deprojected bar radius; col. 6, deprojected bar radius relative to isophotal 
radius; col 7: deprojected near-IR ellipticity of the bar (in parentheses, exp— exponential-profile bar while flat— flat-profile 
bar); col. 8: deprojected bar radius relative to the standard (extinction-corrected) isophotal radius, based on the separated 
bar+disk image; col. 9: deprojected near-IR ellipticity of the bar, based on the same image; col. 10: Fourier near-IR m— 2 
amplitude; col. 11: average inner-limit spiral point angles relative to the bar axis, in degrees. 



Table 4 

Dust-Penetrated Classification 



Galaxy 


Torque Class 


Bar Class 


Spiral Class 


Radius 


Rangc 


m 


1^1 


DP Type 


1 


2 


3 


4 


5 




6 


7 


8 


NGC 908 


3 





3 


24- 


-60 


2 


28.07 


H2/30 














3 


38.66 




NGC 972 


2 


2 


2 


6- 


19 


2 


38.66 


H2 7 2 


NGC 1058 


1 





1 


6- 


12 


2 


25.20 


H2/30 


NGC 1255 


2 


1 


2 


8- 


-27 


2 


31.61 


H2/31 


NGC 1365 


5 


4 


4 






2 




H2 7 4 b 


NGC 1433 


4 


4 


2 






2 




H2a4 b 


NGC 1530 


7 


6 


4 


29 


-58 


2 


33.69 


H2 7 7 


NGC 1808 


2 


2 


1 


18- 


-48 


2 


69.44 


H2 7 2 


NGC 5033 


1 


1 


1 


10- 


-24 


2 


28.07 


H2/31 


NGC 6643 


2 





2 


7- 


-41 






....0 


NGC 6814 


1 


1 


1 


8- 


-36 


4 


34.82 


H4 7 l 














2 


11.31 




NGC 6951 


3 


3 


2 


24- 


60 


2 


45.00 


H2 7 3 


NGC 7217 



















....0 


NGC 7479 


7 


6 


5 


24 


-60 


2 


33.69 


H2 7 6 


NGC 7606 


1 





1 


17- 


-48 


2 


31.61 


H2/30 


NGC 7723 


3 


3 


1 


10- 


-24 






....3 


NGC 7741 


8 


8 


4 


22- 


-34 






....8 



a Col. 1: galaxy name; col. 2: maximum relative gravitational torque class; col. 3: bar torque class; col. 
4: spiral torque class; col. 5: radius range (arcscc) for harmonic and pitch angle fits; col. 6: multiplicity 
of dominant spiral harmonic; col. 7, pitch angle associated with harmonic component m; col. 8: dust- 
penetrated classification (Block & Puerari 1999; Buta & Block 2001). 

b Harmonic and pitch angle classes from Buta & Block (2001). 
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Fig. 1 . Plots of relative Fourier intensity ampli- 
tudes as a function of radius for 12 galaxies where 
bar /spiral separation was needed. Symbols show 
the extrapolations used for our analysis (see text). 
For the stronger bar cases, even terms for m=2 
(solid curve), 4 (dotted curve), and 6 (short dashed 
curve) are shown. For NGC 1255 and NGC 6814, 
only the m = 2 and 4 terms are shown. 



Fig. 2. — Plots of the phase of the m—2 Fourier 
component for the same 12 galaxies as in Figure 
1. The solid vertical lines indicate the radius of 
the bar maximum from Table 2. 
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Fig. 3. — Predicted axisymmctric rotation curves 
for 16 of the sample galaxies, assuming a constant 
mass-to-light ratio at 2.1fim. Each rotation curve 
is normalized to its maximum V m . 



Fig. 4. — Separated bar and spiral images for 12 
galaxies. Three images are shown for each galaxy 
in the following categories: m—0-20 sum (left), 
bar+disk (middle), and spiral+disk (right). Each 
galaxy is identified only in the leftmost of its three 
images. 
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Fig. 5. — Separated bar and spiral force ratio 
maps for 12 galaxies. Three ratios maps are shown 
for each galaxy in the following categories: m=0- 
20 sum (left), bar+disk (middle), and spiral+disk 
(right). Each galaxy is identified only in the left- 
most of its three maps. 



Fig. 7.- Plot of the maximum relative spi- 
ral torque, Q s , versus the maximum relative bar 
torque, Qh- Strong bars and spirals are towards 
the top and right, respectively. 



Fig. 6. — Plots of the average maximum ratio of 
the tangential force to the mean radial force as 
a function of radius normalized to the extinction- 
corrected de Vaucouleurs standard isophotal ra- 
dius, for 16 of the sample galaxies. The sepa- 
rate curves refer to the bar (dashed curve), spiral 
(dotted curve), and the to=0-20 sum image (solid 
curve). 
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Fig. 9. — Plot of the azimuthal angle of the begin- 
ning of the inner spiral, averaged over two or more 
arms, versus the bar strength Qi,. The error bars 
on the angles are average deviations, while those 
on Qb are total mean errors. 



Fig. 8. — Plots of (a) bar-interbar contrast at 
0.7rb or versus Qb] (b) bar ellipticity versus Qb 
(filled circles are eb for the full images, while 
crosses are e b for the separated bar images); 

(c) m—2 Fourier transform amplitude versus Qb; 

(d) bar-interbar contrast versus bar radius rela- 
tive to the (extinction-corrected) face-on standard 
isophotal radius; (e) arm-interarm contrast at 0.25 
times the standard isophotal radius versus Q s ; and 
(f) same as (e), for 0.5 times the standard isopho- 
tal radius. 



Fig. 10. — Inverse Fourier transform contours of 
dominant harmonic superposed on a deprojected 
K s -band image of NGC 1530. 
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